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Tunnel junction between two superconductors is considered in the vicinity of the critical temperature. Super- 
conductive fluctuations above give rise to the noise of the ac Josephson current although the current itself is 
zero in average. As a result of fluctuations, current noise spectrum is peaked at the Josephson frequency, which 
may be considered as precursor of superconductivity in the normal state. Temperature dependence and shape of 
the Josephson current noise resonance line is calculated for various junction configurations. 
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I. INTRODUCTION 

In the vicinity of the critical temperature Tc transport prop- 
erties of metals are strongly affected by superconductive fluc- 
tuations. For example, in the temperature region T - T^, <sc T^, 
where fluctuations are the most pronounced, Drude conductiv- 
ity acquires noticeable Aslamazov-Larkin, Maki-Thompson, 
and density-of-states (DOS) corrections. Many other kinetic 
and thermodynamic coefficients such as magnetic susceptibil- 
ity, heat conductivity. Hall coefficient, and ultrasonic attenua- 
tion are also modified by fluctuations. One may consult recent 
book (Ref.[l]) for exhaustive overview of results and Uterature 
in this field. 

Mostly immediately after the pioneering works on super- 
conductive fluctuations^ai^ it was noticed that analog of the ac 
Josephson effect may survive in the normal state above the 
critical temperature .i^ The latter is also attributed to the for- 
mation of fluctuating Cooper pairs. Indeed, consider weak 
transparency tunnel junction between two superconductors. 
In this case Josephson current is given by Ij{t) - Ic sin(wyf), 
where (jJj = 2eV is the Josephson frequency and current am- 
plitude Ic is proportional to the product of superconductive 
order parameters (^l(R), taken from the left (L) and right (R) 
to the contact area. Above the critical temperature Joseph- 
son current vanishes {Ij) ~ (A^A^) = since order pa- 
rameter is zero in average (A/,(r)) s 0. However, current 
squared {Ij} ~ (A£ArA£Ar), which gives noise of the Joseph- 
son current, is apparently not zero due to nontrivial average 
(A£(s)A/,(s,) of space and time fluctuating order parameters. 
As a consequence, noise power spectrum S j(cl>), defined as 
the Fourier transform of Josephson current-current correlation 
function, shows distinctive peak at the Josephson frequency 
cl>j, which is experimentally an accessible effect. The peak 
height 5 max - S j{a> - wy) is a strong function of T - Tc, 
usually some power law, which makes it possible to detect 
noise signal in the immediate vicinity of the critical temper- 
ature T - Tc ^ Tc- Although this observation was there for 
a long time, the interest to it was recently revived. It was 
stressed^i^ that measurements of the Josephson current noise 
may be especially fruitful in studies of the high-temperature 
superconductivity. Indeed, whether superconductive pairing 
fluctuations exist in the pseudogap regime of the high-Fc ma- 
terials may be probed by the Josephson tunneling. Thus, exis- 
tence of the Josephson effect above Tc may be thought as the 



precursor of superconductivity. 

So far fluctuations of the Josephson current above the 
critical temperature were studied either for the narrow con- 
tacts^i^, taking into account only temporal fluctuations of 
the order parameter, or for the mesoscopic rings^iiS,. We 
find, however, that in the planar geometry of the tunnel 
junction, where spatial variations of the superconductive or- 
der parameter have to be accounted for, peak in the current 
noise spectrum is more pronounced, especially, for the non- 
symmetric junction configurations. Motivated by the ongoing 
experiment&ii and possible applications in probing pseudo- 
gap regime of high-r^. materials, we revisit problem of the 
Josephson current noise above Tc and study noise in the pla- 
nar geometry of a tunnel junction. Within this work we focus 
on the temperature range Gi <sc (T - Tc)/Tc < 1, where Gi is 
the Ginzburg number. In this regime fluctuations can be con- 
sidered as small and can be treated in perturbation theory. The 
natural expansion parameter, which measures strength of the 
superconductive fluctuations, is Gi < I. 

The main results of the present work may be summarized 
as follows: (i) For symmetric wide junctions, when both elec- 
trodes are in the fluctuating regime, and contact area ^ is large 
as compared to the square of the superconductive coherence 
length, y[ > ^l, Josephson current noise spectrum S j(co) has 
a Lorentzian-like shape. The peak height scales in tempera- 
ture as 5 max [x-t] depends quadratic ally on both tun- 
nel conductance of the junction gj and the Ginzburg number 
Gi. For the lowest temperature T — Tc - Gi Tc, which is al- 
lowed by the applicability of the perturbation theory, strength 
of the noise is given by 5max = {nl6A){g].Tcle^){^llJ[). Of 
course, experimentally, noise is maximal right at the tran- 
sition T - Tc', however, in this case it is very difficult to 
make any quantitative predictions theoretically. Thus, 5 max 
gives an order of magnitude estimate, (ii) For the narrow, 
~ symmetric junctions we find also a Lorentzian-like 
shape of S jiai), which is again quadratic in both gr and Gi; 
however, temperature dependence of the peak height is dif- 
ferent 5 max <^ estimate for the noise power at the 
most vicinity of the transition is 5 max - {Gi/^n)(g^Tc/e-^). 
(iii) In the case of non-symmetric junctions, when one elec- 
trode is akeady superconducting while another is fluctuating, 
noise has Lorentzian form. The temperature dependence for 
the peak height in this case is the same as for wide symmetric 
junction, which, however, appears already in the first order of 
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the Ginzburg number and contains large pref actor ln^(As jTc) 
(where is the superconductive gap), (iv) Corrections to 
the current noise above Tc are not exhausted by the Josephson 
current contribution only. In addition, superconductive fluctu- 
ations deplete normal-metal DOS at the Fermi energy, which 
changes tunnel conductance. The latter translates into the cur- 
rent noise correction 5dos(<j^) via fluctuation dissipation the- 
orem (FDT). This eff'ect is linear in gj and Gi, logarithmic 
in temperature ^dos ^ ^'^(r^)' ™ opposite sign as 

compared to the Josephson current contribution. 

The rest is organized as follows: in the next section (Sec.HIll 
we present in a concise form our technical method, Keldysh 
nonlinear cr-model, which will be used through out the paper 
in calculation of the current noise power. This formalism was 
elaborated in Refs. [13 and [TH and found to be very useful 
and powerful in many applications. In the Sec. |lll]we calcu- 
late density-of-states and Josephson current contributions to 
the noise spectrum above T^.. The results of the work together 
with further discussions are summarized in the Sec.|IV] Num- 
ber of technical points are delegated to the Appendixes IV Al - 
IVCl 



II. FORMALISM 

Consider voltage biased tunnel junction of two supercon- 
ductors above the critical temperature. Within cr-model for- 
malism tunneling between L and R reservoirs of a junction is 
described by the action. 
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where gj is the junction tunnel conductance and Q,l{R) are the 
Green's functions describing electron system in the electrodes 
(hereafter - kg - 1). Both Q,l{R) are 4x4 matrices in the four 
dimensional Keldysh i8>Nambu space. Matrix H = o-q i8> t-, 
where cr,, t, for / - Q,x,y,z, are the sets of Pauli matrices 
acting in the Keldysh and Nambu subspaces correspondingly, 
and symbol ® stands for the direct product. Matrix V is the 
source term having standard structure in the Keldysh space. 
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Diagonal elements of V are directly related to the classically 
applied voltage V '(f) = ^Vt, while V^it) is just its quantum 
component. This terminology stems from the Keldysh contour 
- terms classical and quantum imply the symmetric and anti- 
symmetric Unear combinations of the field components resid- 
ing on the forward and backward parts of the Keldysh contour, 
respectively.''* Finally, trace operation Tr{. . .) in Eq. ([TJ as- 
sumes summation over the matrix structure as well as time and 
spatial integrations. The origin of phase factors exp[+;Hy] 
in Eq. ([Til is from gauge transformation, which moves differ- 
ent electrochemical potentials of electrons in the leads from 
the Green's functions to the tunneling term. Dynamics of the 



Green's functions is governed by the cr-model action, 
iSAQL,QR]^-Y, ^TrjAA} 
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where Va is the bare normal metal density of states at the Fermi 
energy. Da is the diffusion coefficient, /!„ is the superconduc- 
tive coupling constant, and S = <Tx®tq. The matrix supercon- 
ductive order parameter Aair, f) is 
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Action ([3]) is subject to the nonlinear constraint Ql - I - Phys- 
ical quantities of interest are obtained from the action via its 
functional differentiation with respect to the appropriate quan- 
tum source. For example, tunnel current is found from the 
equation 
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Z[V] = j D[eje'^ieL,G«]^ (5) 



where S[2z,, G«] = S0-+S7. Corresponding noise power spec- 
trum is defined as 
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The procedure of extracting physical observables, outlined 
above, is rather general within Keldysh technique. However, 
for the problem at hand, information encoded in the actions 
([TJ and ([3]) is excessive. Indeed, S^,- describes not only dy- 
namics of the order parameter A but also contains explicitly 
electronic degrees of freedom in the form of the Q matrices, 
which complicates further analysis. Simplification is possible 
realizing that dynamics of Q is fast as compared to that of A. 
The latter is governed by the time scale tq ~ l/T, while the 
former by ta ~ l/(T - T^), and noticeably ta » tq when 
T ~ Tc- Under this condition, one may integrate out fast elec- 
tronic degrees of freedom from action (O and find an effec- 
tive theory, which describes space and time fluctuations of the 
superconductive order parameter only. This program was re- 
alized for Eq. (O in the recent worbi^ and we will follow here 
the same route in dealing with the tunnel term S7[y]. 

Let us outline essential elements of the method. Having 
interest in the effects of superconductive fluctuations, it is rea- 
sonable to start from the normal metal state with the Green's 
functions Qur-, = Qn given by 
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which minimizes action (O for A^ = 0. One treats then A^ 
in perturbation theory on top of Qff. Technically this program 
is realized in several steps. At the first stage one projects Q- 
matrices as 
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where W"^,(r) carries information about fast electronic de- 
grees of freedom. Matrix W is parametrized by the two com- 
plex fields c^^i{q) and Cce'iq) — Cooper modes, which will be 
integrated out eventually. It is convenient to choose 



W" ^Ro'W" oR- 



(9) 



with 



= T+ + c::,T_) ® cr+ + (c:,,t+ + t_) ® (r_, (lO) 
where t± - (tj. ± iTy)/2, cr-t = (ctq + cr-)/2 and 
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One brings then Eq. ([8]l into action ([3]) and expands So-[2a] 
^a-[W",Aa] to the second order in the Cooper modes W" = 
' ^L' 1 (details of this procedure are provided in the Ap- 
pendix |VA]i. One finds then that to the leading order in the 
coupling TrjQA), Cooper modes are connected to the super- 
conductive order parameter according to the relations 

where we have introduced retarded(advanced) Cooperon 
propagator, 

1 



(13) 



Daq^ + i{€ + e'y 
and the form factors, 

\'^^,(q) = -2[Ai^,(q) + F,Al^,(q)], (14) 

Knowing relations (fTSl i Gaussian integration over the Cooper 
modes is straightforward, 



D[W"] sxpiiSAW, A„]) = exp(iSefF[A]). (15) 



The corresponding quadratic form Str[ W, A^] should be taken 
from Eq. (l47T i and one finds as a result, 

SefF[A] = 2 2v«Tr{A^L^iA«}, A,^ = (A^',A^). (16) 



The propagator £ '(q', cj) governs superconductive order pa- 
rameter dynamics. It has typical bosonic structure in the 
Keldysh space 



(17) 



with 



8r, 



-{Daq^ + Tcl + iaj), (18) 



Zj^(q,a)) = B^[Si^^{q,a))- iij^^{q,(jj)\, 

and tgl — ^/^{T - Tea) and B^^ = coth(w/2r). 

Noticeably, effective action (fT6] l is much simpler than the 
original one [Eq. Q]. However, what is important to empha- 
size, is that Seif captures correctly all the relevant low energy 
excitations of A„{r, t). After these technical preliminaries we 
turn now to the applications of the general formalism based 
on the effective action SefF[A]. 



III. CURRENT NOISE ABOVE T, 

A. Tlinnel current noise 

The first apparent effect of superconductive fluctuations is 
modification of the normal metal density of states. Being flat 
in the normal state, v(e) acquires strong energy dependence 
in the vicinity of with a dip around Fermi energy^i^ The 
latter suppresses tunnel conductance of the junction, which 
influences tunnel current and as the result its noise. Supercon- 
ductive fluctuations correction to the tunnel current was stud- 
ied in Ref. [ItI Here we calculate corresponding correction 
to the noise. Although the result of this calculation follows 
immediately from the fluctuation-dissipation relation it is still 
useful to see how it appears within the cr-model approach. To 
this end, assume non-symmetric tunnel junction: let us say 
that left electrode is in its normal state, while the right one 
is in the fluctuating regime. To calculate noise power, one 
uses general definition [Eq. Q] and inserts Qi - Qn and 
Qr ~ Qn[1 + iW - (Ref.[ll) into the tunneling part 

of action ([TJ. After the difl'erentiation, which is done with the 
help of the formula 
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= ±i6(t - t')fexp[+ieVtS] , (19) 
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where Y - <t_^ ® r., one finds for the noise 



5(w) = Ss(co) + 5dos(w). 



(20) 



Here 



Ss{co) = 2gTTY,^coth^, (21) 

± 

with u± - eV + CD, is just the Schottky formula for the noise 
in the normal tunnel junction, while the corresponding fluctu- 
ations correction is 

d{t-t')[SAt,t') + s.it,t')y"'^'-''\ 

CO 

(22) 

where 

S^{t,t') = Tx[Qf,{e)RA{'W,AqW,.,{-q))) 

^fteA,(e")te^'''^""'''^e'='''"'"^*'"''^} . (23) 

Quantum averaging in Eq. ( [23] ), denoted by the angular brack- 
ets ((. . .)), should be performed with effective action ( [T6l ), 
namely, {{. . .)) = /D[A] . . . exp (/Seff [A]). Recall that fluc- 
tuation matrix 'W is expressed through the Cooper modes c^f- 
and Cff', which are functionally dependent on the order pa- 
rameter A via Eq. ( [T2b . The notation 5dos in Eq- ( l22b and its 
actual relation to the density-of-states suppression are moti- 
vated in Appendix [VB] The linear in 'W term in Eq. ( [23] ) is 
not written explicitly since it does not contribute to the final 
result. The final comment in order of Eq. ( l22l i is that traces 
of iS± functions allow rather simple and convenient diagram- 
matic representation shown in Fig.[T^. 
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At this point one calculates the product of 'W matrices in 
Eq. ( |23] l and performs Gaussian functional integration over 
the fluctuating order parameter using Eqs. (fT2l i and (fT6] l. The 
resulting averages are 



^Kiq, e-e') + F^^LRiq, e - e') + F^LAjq, e - e') 

{(c^AqK'A-q)}} = (2//V) 
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Next few steps are conceptually simple, (i) One traces 
Eq. ( |23T l over its matrix structure first and then performs 
time Fourier transforms in Eq. ( l22l i J d{t - f')e'(f-f"±«*'K'-'') - 
2nS{e — e" + eV), which removes e" integration, (ii) Ob- 
serve that for the e' integration, term containing F^L/Sji, e - 
e') in the average {{cc*}} and term containing F^£jR{q,£ - 
e') in the average {{cc*}} do not contribute to S± as be- 
ing integrals of purely advanced and retarded functions, re- 
spectively. As the result, one takes {{cc*}} + {{cc*}} = 
2/Im{(cc*)). Finally one changes momentum sum into the 
integral 2^ — > / d~q/4-n^, assuming that the electrodes are 
quasi-two-dimensional films, and introduces dimensionless 
variables x = Dq^/T, {/ = (e - e')/T, and z ^ (e + e')/4T. 
After these steps Eq. (l22l i becomes 
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Here rj - \ /TcTcl, and we introduced Ginzburg number Gi - 
1 /vD. After the remaining integrations (see Appendix lV Cl for 
details) one finds as a result, 

Sdos(oj) = -^grrinj^^j (27) 

X Vcoth^Im^t'lfi-i^ 



where (/(''^(z) is the first-order derivative of the digamma func- 
tion. Close look on Eq. dZTl l allows us to rewrite it in the form 



Sdos((^) = « ^ feos(M±) coth ^ , 



(28) 



where /dos is the tunnel current correction calculated in 
Ref.[l3, which is a priori expected result from FDT. 

In complete analogy one can calculate corresponding cor- 
rection to the noise for the symmetric junction when both elec- 
trodes are in the fluctuating regime. In this case Green's func- 
tion matrix Qi has to be expanded in fluctuations lV also and 
one faces diagram shown in Fig. [TJi. The result of the calcu- 
lation can again be cast in the form of Eq. (|28T l. where /qos 
should be replaced by the appro pria te second-order fluctua- 
tion correction known from Ref. [l3 Furthermore, if one is 
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FIG. 1: Superconductive fluctuation contributions to the current 
noise. Diagrams (a) and (b) correspond to the effects coming from 
the fluctuations in the density of states for the non-symmetric and 
symmetric junctions. Diagrams (c) and (d) are the fluctuating Joseph- 
son current contribution for the symmetric and superconductor- 
fluctuating metal junctions correspondingly. Ladders represent 
Cooperons, Eq. l il3t , wavy lines stand for the fluctuations propaga- 
tor, Eq. J18t . and crossed boxes depict tunnel conductance gj. 



able to calculate /qos completely, meaning to all orders of 
perturbation theory, then for the noise of the tunnel current 
Eq. ( |28] ) can be considered as the exact result, which is again 
consequence of FDT. 



B. Josephson current noise 

Contribution to the noise spectrum coming from the 
Josephson effects is very much different than that of density 
of states. First of all there is no simple FDT relation similar 
to Eq. ( l28T l. Secondly, the physical mechanism, which leads 
to the noise, is different. Probably the simplest way to see this 
is to start from the definition of the current in Eq. As- 
suming symmetric junction configuration, one expands then 
each Green's-function matrix to the linear order in fluctua- 
tions Ql(R) — > iQNWiiR) in the tunnel part of action which 
gives for the current. 
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To proceed further, we will simplify Eq. ( |29] l, exploring sepa- 
ration of the time scales between electronic and order param- 
eter degrees of freedom. Indeed, one should notice that as it 
follows from Eq. ( fTSb relevant energies and momenta for the 
order-parameter variations are Dq^ - oj ~ r^J^, while the rel- 
evant fermionic energies entering the Cooperon in Eq. ( fT3l ) 
are e ~ e' ~ \jT . As a result, nonlocal relations between 
Cooper modes and order parameter in Eqs. ( fTOl i and ( fTSI ) can 
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be approximated ai 
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where 0f are the retarded (advanced) step functions. Phys- 
ically Eq. ( [30l l implies that Cooperon is short-ranged, having 
characteristic length scale - ^IDITc, as compared to the 
long-ranged fluctuations of the order parameter, which prop- 
agates to the distances of the order of ^gl = VDtgl s> ^q. 
Thus, relations (fT2l i are effectively local, which simplifies fur- 
ther analysis considerably. Equations (l30l l allow us to trace 
Keldysh subspace in Eq. ( [29] l explicitly to arrive at 
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where we have used Eq. ( fT9] l and wrote trace in the real space 
representation (note that Tr{. . .) here does not imply time f 
integration). Changing integration variables t[ - t - and 
f3 - t - T], and rescaling 77, /i in the units of temperature Trj 
T], Tji — > ju, one finds forEq. dSTT i an equivalent representation. 
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where we used equilibrium fermionic distribution function in 
the time domain F, = -;T/ sinh(7rrf). The most significant 
contribution to the above integrals comes from rj ~ ji < \. 
At this range ratios {r],fi}/T change on the scale of inverse 
temperature, while as we already discussed, order-parameter 
variations are set by f ~ tql '» l/T. Thus, performing 77 and fi 
integrations one may neglect {rj, iJ.]/T dependence of the order 
parameters. As the result we find 



J ^[A^'(r,OAr'(r,Oe-'"" 
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Finally we are ready to calculate corresponding contribution 
to the current noise. One brings two currents from Eq. ( |33] | 
into Eq. (|6]l and pairs fluctuating order parameters using cor- 
relation function. 



«A^'(r, t)Af(r', f')»A = ^^ab^dr -r',t- t') , (34) 



which follows from Eqs. ( fTSI l and (fTST l. As a result, Josephson 
current correction to the noise of wide symmetric junction is 

(35) 

where u-t - u + toj. Corresponding diagrammatic represen- 
tation of Eq. ( I35] ) is shown in Fig.[T];. Remaining integrations 



in Eq. (|35] | can be done in the closed form (see Appendix lV CI 
for details), providing 



S j(oj) = > — A'(w±tgl) 



N{z) = -T In Vl + zV4. 
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Analogous calculation in the case of the narrow symmet- 
ric junction, which is obtained from Eq. ( [35] ) by replacing 
i^Kif, t) — > ^k{Q, t) and removing spatial integration, gives 
for the noise spectrum (see details in Appendix lV Cb 
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In a similar fashion one may consider non-symmetric tun- 
nel junction. Assume that one of the electrodes is in the deep 
superconducting state, with well defined gap in the excitation 
spectrum A5 , while the other is in the fluctuating regime. We 
set then one of the Qa matrices to be superconductive Green's 
function Ql — Qs^ where 



Qs 

F - Ft., and 
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while expanding the other one in Cooper modes Qr 
IQn^r. The resulting expression for the current reads 

'^^'^ = -^^Tr{.'-^^&.-'-^^GMW«} . (40) 

Following the same steps as in the case of the symmetric junc- 
tion, carrying out differentiation with the help of Eq. ( fT9l ) and 
tracing consequently Keldysh and Nambu subspaces and per- 
forming time integrals, one finds for the current. 
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where we assumed that As » Tc- Squaring Eq. (HTT i and av- 
eraging over the order-parameter fluctuations with the help of 
Eq. (|34l i, we get 

(42) 

Performing the remaining integrations, one finds noise spec- 
trum of the non-symmetric junction [see corresponding dia- 
gram in Fig.[T]J], 
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FIG. 2: (Color online) Shape of the Josephson current noise spectral 
lines in the vicinity of the resonances z = ai±TGL- 



Spectral line shapes for Eqs. ( [36] l. dJTj i and ( |43] l are plotted in 
the Fig. H 



IV. DISCUSSIONS 

We have considered effects of superconductive fluctuations 
on the current noise in tunnel junctions above the critical tem- 
perature. Several contributions were identified. The sim- 
plest one originates from the fluctuation suppression of the 
density of states. This effect gives negative contribution to 
the current noise, which is only logarithmic in temperature 
5 DOS 111(2" - Tc), whereas dip in the density of states at 
the Fermi energy has much stronger temperature dependence 
6v{Q) cc {T - TcY^- Somehow current and its noise get sup- 
pressed weaker than the density of states itself. Another in- 
teresting point is that current noise is strongly modified only 
at the characteristic voltages eV ~ T^, while corresponding 
feature in the density of states appears at energies e ~ T - T^, 
see Eq. ( |5T] i. It turns out that higher-order fluctuation effects, 
similar to that shown in Fig. [T}5, restore additional structure 
of the noise signal at eV ~ T - T^. Correction ^dos is lin- 
ear in Gi and in tunnel conductance gr- This is in contrast to 
the Josephson current contribution to the noise. The latter is 
quadratic in fluctuations and in tunneling, and enhances noise 
at the frequencies in the vicinity of the Josephson frequency 
cl>j. The peak at w = is well defined and is strongly tem- 
perature dependent, which makes it possible to detect it exper- 
imentally. We have found that depending on the junction con- 
figuration: symmetric or non-symmetric and narrow or wide, 
noise resonance line has different shapes in the frequency do- 
main Fig.|2]and different temperature dependencies. 

Closing this section we should mention that in the field 
of fluctuating superconductivity one usually identifies three 
types of fluctuation corrections. Apart from density of 
states, there are also so called Aslamazov-Larkin (AL) and 
Maki-Thompson (MT) terms mentioned in Sec. I. It is quite 



natural to ask how AL and MT processes modify current 
noise and how they can be identified within cr-model for- 
malism. As an attempt to answer, one should recall that 
in addition to the simple tunneling term Sr[y], consid- 
ered in this work, one may have yet another one /S/i[y] = 
{7TgA/l6e^)TT{{e'^^ Qie^'^^ Qk)"}, which was neglected. It 
corresponds to Andreev processes, and gA is Andreev con- 
ductance. Using Sa[V], instead of S7-[y], one may follow the 
same routing expanding Q-matrices in fluctuations W to ob- 
tain additional contribution to the noise. However, among all 
the terms emerging in perturbative expansion, separation on 
AL and MT contributions becomes ambiguous. Nevertheless, 
the problem is very interesting and requires further studies. 
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V. APPENDICES 

A. Fluctuations expansion 

Within this section we show in details how the transfor- 
mation from Eq. Q to Eq. ( fT6] l occurs. We start fluctua- 
tions expansion by taking Q » 2^(1 + " W~/2) and 
bringing it into the So- (here and in what follows subscript 
a in Qa matrix and all other elements will be suppressed 
for brevity). For the trace of the gradient term we find, 
Tr{(V02) = -TtIQn'^WQn'^W] = -Trl'lVV^'VV), where we 
employed anti-commutativity relation [Qm, W]+ = and non- 
linear constraint - 1 . Using an explicit form of the 
matrix [Eq. ( fTOl il and tracing the product of two "W over the 
Keldysh i8> Nambu space, we obtain 

Tr{(V02} = 2Tr{Dq\c:Aq)Ce'e(-q) + c:Aq)Ce'e(-q)]} ■ 

(44) 

The time derivative term in the So- produces contribution 
Tr{&d,Q} = ^Tr{eR,M,Any,,,'W,,,}, while Unear in W 

part traces out to zero [here we used Qn{£) = R^ARf^ with 
A = (T- ® Tj and substituted d, — > -ie]. Observing that 
RgERirA - CT; (g) To, one finds 

Tr{B5,e} = ^Ti{(e+e')[cUqM-q)-cUqrc.'.(-q)]]- 

(45) 

For the coupling term between Cooperons and A, to the lead- 
ing order, we have Tr{Ag} = Tr{ReA,^e,Re,A'W,,'l which 
translates into 

Tr[AQ] = -m{[AtAq) + F,Al_Aq)]cl,,{-q) (46) 
+ [Af,{q) + F,Al*_Aq)]c,,,{-q) - [AtAq)' 
F,,AlAq)K,,(-q) - IK'-Aq) - F,,Al%(q)]c,A-q)} ■ 
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Combining now Eqs. (I44li-(|46]| all together and bringing them 
back into the Eq. (O, we wind for the quadratic in Cooperons 
part of action SAW, A] = Sj;[c, A] + S^[c, A], where contri- 



butions from the retarded c and advanced c Cooperons read 

as 



/S^[c, A] = -—Tr{c:ADq' - i(e + e')]c,,, + 2[A^ 



-£ t 6-6' 



(47a) 



(47b) 



r 



At this stage we are ready to perform integration over the 
Cooperon modes. Assuming that configuration of the order- 
parameter field is given, one varies Eq. ( |47] ) with respect to c* 
and c*, and obtains stationary point equations 5SJ^/5c*^, = 
and SEi^j./6c*^, - 0. The latter are easily solved by Eq. (fT2l i. 
Since the value of the Gaussian integral is equal to that taken 
at the saddle point, one brings Eq. ( fT2] i into Eq. ( |47] | and after 
some straightforward algebra finds Eq. (fTSI l. Further details 
can be found in Ref. [isL 



B. Relation between 5dos and <5y(e) 

The purpose of this section is to demonstrate explicitly that 
S DOS indeed originates from the DOS eff'ects, which was hid- 
den in the technical details of Sec. To this end we cal- 
culate temperature dependence of the 6v(e) within Keldysh 
technique. This illustration is useful for the sake of compari- 
son with the known results obtained previously from the tem- 
perature Matsubara technique^. 

Within cr-model energy dependent density of states is ex- 
pressed in terms of Q matrix in the following way: 



v(e) = -TilQ^Qeel 



(48) 



Setting Q - Qf^ one recovers bare normal-metal density of 
states v(e) = v. To account for the fluctuations on top of 
the metallic state, one expands Q in Cooper modes W to the 
quadratic order and averages over A fluctuations with the ef- 
fective action from Eq. ( fT6b : 



6v{e) 



-^Tr[«c, 



(49) 

Observe that this is precisely the same combination of the 
Cooperons, which enters ^dos in the Eq. ( l22l i. thus they 
have common origin. Furthermore, it is easy to show that 



6+IU/2T 



- F 



e-«±/2rJ 



from Eq. (l24b , density-of-states correction becomes 



Using averages 



6v(e) 



In {D(f- - lie + iojf 



(50) 



where we set e' - e - u. Here one meets the convenience of 
the Keldysh technique, which allows us to get physical quan- 
tities avoiding analytic continuation procedure. Using explicit 



form of fluctuations propagators from Eq. (fTSI ) and perform- 
ing frequency and momentum integrations, one finds in the 
quasi-two-dimensional case. 



5v{e) 



Gi_ 
'16 



J - T, 

where dimensionless function is 

dx 



T{z) = Re 



Jo 



(1 +x){\+2x-2izy 



(51a) 



(51b) 



In agreement with Ref. 111611 dip at the Fermi energy is 5v(0) oc 
{T - TcY^, while at large energies ercL s> 1 density-of-states 
correction recovers its normal value according to 6v{e) oc 
(r,/e)2ln(eTGL). 



C. Integrals for 5 DOS (w) and 5 J 



(I) Transformation from Eq. 
culation of the integral. 



to Eq. ( |27] i requires cal- 



Jo Jj-(X 



dydzRc- 



F- 



+i<±/2T 



z-iuIlT 



' [{x + 77)2 + y^{x + iy - Aiz)^ ' 

(52) 

One performs y integration first. 



- n I dx I 

Jo J-<x 



dzRe 



'z+«i/2r ■ 



-IU/2T 



{x + t])(t] + 2x- 4iz)^ ' 



(53) 



Since 77 <«c 1 and relevant z ~ 1 one may safely approximate 
{rj + 2x - 4iz) ~ 2{x - 2iz). Then expanding F~ into the series 
= 2 J]„z/(z- + zl), with z„ = 7i(n + 1/2), interchanging 
order of summation and integration and recalling definition 
of the «th-order derivative of the digamma function i/''"l(z) = 
(-!)"+'«! 2,7=0 1 /(« + ^)"^'' one finds that 



f 

•J — cx 



dz 



'z±«±/27- 

(x - 2/z)2 



I 



1 



iu± 



— ± — + — 

27T^ \2 27tT 271 



(54) 



Remaining x integration can be taken with logarithmic accu- 
racy, ignoring x dependence of the digamma function since 
only X < 1 contribute significantly, which eventually gives 
- In rj. Combining all together, one finds 



(55) 
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which in combination with Eq. ( |26] ) resuhs in Eq. dZTl l. 

(II) Transition from Eq. ( [35T l to Eq. ( [36] l is performed in the 
following way. As the first step one finds Keldysh component 
of the fluctuation propagator in the mixed momentum/time 
representation LKiq, - J ^Kiq, to)e^"^'dto/2n, which gives 

LKiqj)^-—^ , 7Cg = {^CLqf + l- (56) 

One inserts then ^Kir, t) - J J^xiq, t)e"'''dq^ /4k into Eq. dSST l. 
integrates over r, introduces dimensionless time r = I/tgl, 
and changes from q to x integration dq^ - dx/^^^, which 
gives all together 



which defines N(z) function in Eq. 

(Ill) Calculation of Eq. ( l37l i is completely analogous. 
Noticing that £;(-(0, f) = J iLjf(q', f)c/(7^/47r and transforming 
to the dimensionless units r - T/tgl and x = (^gl^)^ + 1, we 
have 



± 

dT I 



dxdx' 



KX' 



-(z+x')|T|-iZ±T 



(59) 



T -T, 

+ 00 r^+QO 



X 



Ji x~ 



where z± — co±TGh- After t integration one is left with 

4dx 

li 



r 



x(4x^ + zi) ' 



(58) 



After T integration one is left with 



rp°° t/;>it/;>i' ?i + ;>i' _ 2 c/;^ ^^^^^^ f ^ 

JJl XX' (k + x')^ +z\ z± Ji X \ z± ) 

(60) 

which after the integration by parts reduces to 2M{z), with 
M{z) function defined by Eq. ( [37] i. 
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